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5-H , Abstract 

We study the classical geodesic motions of nonzero rest mass test particles and photons in (3 + 1 + n)- 
dimensional warped product spaces. An important feature of these spaces is that they allow a natural 
decoupling between the motions in the (3+ l)-dimensional spacetime and those in the extra n dimensions. 
Using this decoupling and employing phase space analysis we investigate the conditions for confinement 
of particles and photons to the (3 + 1)- spacetime submanifold. In addition to providing information 
] regarding the motion of photons, we also show that these motions are not constrained by the value of 

the extrinsic curvature. We obtain the general conditions for the confinement of geodesies in the case of 
pseudo-Riemannian manifolds as well as establishing the conditions for the stability of such confinement. 
These results also generalise a recent result of the authors concerning the embeddings of hypersurfaces 
with codimension one. 



1 Introduction 



An intriguing idea in modern cosmology is the possibility that the universe may have a higher number of 
' dimensions than the classically observed (3+1). There has been a number of motivations for this idea, mainly 
related to attempts at constructing a fundamental theory of physical interactions. These range from the 
original attempts at unification of gravity with electromagnetism by Kaluza and Klein [1] to String/M-theory 
[2,4, 5]. ^ 

An immediate task within this higher dimensional framework is how to explain the four dimensionality 
of the observed universe. An original idea - which dates back to the work of Kaluza and Klein - assumes 
these extra dimensions to be compact and very small. In most recent braneworld models, on the other hand, 
stringy effects are invoked to argue that in low energy regimes particles are restricted to a special 3 + 1 
brane hypersurface, which is embedded in a higher-dimensional bulk, while the gravitational field is free to 
propagate in the bulk [6] . 

The intense recent interest in the higher-dimensional scenarios has also provided strong motivation to 
look for geometrical mechanisms which could allow confinement and this has led to the investigation of the 
so-called warped product spaces [7] and their geometrical properties. 

Originally most braneworld and higher dimensional studies concentrated on (4 + 1) scenarios, i.e. co- 
dimensional one models. However, given the prediction of the string theory, according to which spacetime 
is 10-dimensional, the study of higher co-dimensional models is urgently called for. Increasing effort has 
recently gone into the study of such models (see for e.g. [9]). 

An important question concerning such models, specially from an observational point of view, is the 
behaviour of the geodesies in these models, and in particular the relation between the geodesies of the 
higher-dimensional space and those belonging to the hypersurface. A great deal of effort has recently gone 
into the study of geodesic motions in five-dimensional spaces [12]. Recently the authors studied the ability 
of five-dimensional warped product spaces to provide a mechanism for geodesic confinement on co-dimension 



one hypersurfaces, purely classically and based on gravitational effects [11]. The aim of this article is to 

generalise that analysis to the case of co-dimension n warped product spaces [10]. Wc do this by studying 
the classical geodesic motions of nonzero rest mass test particles and photons in (3 + 1 + ri)-dimensional 
warped product spaces. We show that it is possible to obtain a general picture of these motions, using 
the natural decoupling that occurs in such spaces between the motions in the n extra dimensions and the 
motion in the AD submanifolds. This splitting allows the use of phase space analysis in order to investigate 
the possibility of confinement^ and the stability of motion of particles and photons to submanifolds in such 
(3 + 1 + ri)-dimensional spaces. 

The paper is organised as follows. In Section 2 we write the geodesic equations for warped product spaces 
and consider the parts duo to 4D and the higher dimensions separately. We then show that the equations 
that describe the motion in higher dimensions decouples from the rest. We proceed in Section 3 to find 
the mathematical conditions that must be satisfied by the warping function / in order for timelike and 
null geodesies of the higher-dimensional space to be confined to M^. In Section 4 we rewrite the geodesic 
equations in the higher dimensions as an autonomous 2n-dimcnsional dynamical system. We then employ 
phase plane analysis to study the motion of particles with nonzero rest mass (timelike geodesies) and photons 
(null geodesies) respectively. Such qualitative analysis allows a number of conclusions to be drawn about the 
possible existence of confined motions and their nature in the neighbourhood of hypersurfaces. In Section 
5 we give an analysis of the motion in the extra dimensions by reducing the problem to the motion of a 
particle subjected to the action of an effective potential and illustrate the method with an example from the 
literature. We conclude in Section 6 with some final remarks. 



2 Warped product spaces and geodesic motion 

In general a warped product space is defined in the following way. Let {M"^,h) and {N",k) be two Rie- 
mannian (or pseudo-Riemannian) manifolds of dimension m and n, with metrics h and k, respectively. Let 
/ : — > R be a smooth function (which we shall refer to as the warping function). We can then construct 
a new warped product Riemannian (pseudo-Riemannian) space by setting M = M™ x A^" and defining 
a metric g = e^-^h ® fc. In this paper we shall take to = 4 and identify Af^ with the (3-1-1 )-dimensional 

spacetime, a four-dimensional Lorentzian manifold with signature (H ). Therefore the class of warped 

geometries which we shall consider can be characterised by the following line element 

dS'^ = e'^fhafsdx^'dx'^ - kabdy^dy^, (1) 

where / = f{y^, ■•.,?/"), = ha^ix) and kab is in general an n-dimensional Riemannian metric^. 
Let us now consider the equations of geodesies in the {n + 4)-dimensional space M 

dA2 dX ' dA 

where A is an affine parameter and F^^ are the Christoffel symbols of the second kind defined by I'^c ~ 
^9^^{9db,c + 9dc,b — gBCo)- After some algebra it is not difficult to show that the geodesic equations 
for the 4D part can be written in the form 

rfA2 + d\ dX ' 

where 

dy" dy^ rfa;" dy'' 

^ dX dX dX dX 

and ^''^r[^^ = 1.9'^'' {9va,i3 + 9vp,a — 9a0,v)- Likewise the extra-dimensional part of the geodesic equations 
(2) may be written as 

d^y" , (n) dy' dy 



(4) 



dA2 '"=dX dX 

where 

_ p„ dx- dx^ dx- dy'' 



'"'nc'^^ = r (5) 



'^Throughout by confinement of photons we mean to say that their motion is constrained to he in the brane, rather than 
being bounded. 

^Throughout capital Latin indices take value in the range (0, 1, ...(3 + n)), lower case Latin indices take values in the range 
(4, ...(3 -|- n)) while Greek indices run over (0, 1, 2, 3). Thus, the coordinates of a generic point P of the manifold M will be 
denoted by = (s™, where x" denotes the 4D spacetime coordinates and j/" refers to the n extra coordinates of P. 



We now assume that our spacetime corresponds in this scenario to a particular submanifold defined 
by the n equations = =constant. The geometry of is then determined by the induced metric 

ds^ = gaj3{x,yl,...,y^)dx°'dx^ . 

Therefore the quantities ^^^T^^ which appear on the left-hand side of Eq.(3) may be identified with the 

ChristofFcl symbols associated with the metric induced on the leaves of the foliation defined above. 
For the warped product space (1) the quantities (j)^ and i/'" reduce respectively to 

r = -2f,ax^y\ 

and 

where /.a = a dot denotes differentiation with respect to the affine parameter A and we are using the 

notation />" = k°''^ffi. Note that if the warping function / is constant the equations (2) which describe the 
geodesies of the entire warped product M separates into the geodesic equations for the two submanifolds 
and . On the other hand if we restrict ourselves to the class of timelike or null geodesies of M we 
have the following first integral 

e-'^K^x'^xl^ - kabry'' = e, (6) 

where e = 1 in the case of particles with nonzero rest mass and e = in the case of photons. With the 

help of the above first integral it is now possible to decouple the motion in the extra dimensions from the 
spacetime. Thus the geodesic equations (2) finally take the form 

x^ +(4) r'^^x'^xf' = -2lax^r, (7) 

r+ncy''r = -rie+kbcy'r)- (s) 

These equations constitute a system of second-order ordinary differential equations which can, in principle, 
be solved once the function / = f{y^, y") is given. 



3 Confinement of the motion in four-dimensional spacetime 

In this section we wish to investigate the possibility of confinement of massive particles and photons in the 
spacetime manifold Af*. This amounts to finding the mathematical conditions required to be satisfied by 
the warping function / such that timelike and null geodesies of M coincide with those corresponding to M^. 
In general such submanifolds of a Riemannian (or pseudo-Riemannian) manifold M are referred to as totally 
geodesic. In the case where M is a Riemannian space, then according to a theorem of differential geometry 
the submanifold is totally geodesic if and only if all normal curvatures (or extrinsic curvatures) of 
vanish [13]. If the geometry of M is pseudo-Riemannian (in our case M is a Lorentzian warped product 
space) then, as we shall show below, the theorem is still valid, however the confinement of null geodesies 
does not depend on the value of the normal curvatures of the submanifold M^. From a physical standpoint 
it is interesting that for the large class of warped product spaces defined by the equation (1) the motion of 
photons is not constrained by the extrinsic (or normal) curvature of the spacetime. This result generalises 
a recent result of the authors concerning the embeddings of hypersurfaces with codimension one [11]. 

To show this, consider equations (7) and (8). Let 7 be a timelike (or spacelike) geodesic curve of the 
submanifold M"^. Since 7 € we must have = 0, which implies that (7) is identically satisfied. Now 
if /'"(j/o, ...y" ) = 0, then (8) also holds. Therefore Eq. (2) is satisfied, implying that 7 is a geodesic of M. 
Conversely, if any geodesic 7 of M with parametric equations {x" = a;"(A),y° = y^) is a geodesic of 
then from (8) we must have f'°'{yl, ■■■yo) = 0. Let us now calculate the normal curvatures of M^. From (1) 
we see that the vectors A''(a) = (where no summation is implied over the index a) are normal to the 

submanifold M^. Let 7 be a curve of with tangent vector given hy V = (^^, 0). The normal curvature 
of M^, at a point p G M^, in the direction of N(^a) is given by the inner product fla = g{N(^a), "^) at p, 
where denotes the covariant derivative of V with respect to the Levi-Civita connection determined by 
g. After computing from (1) we can easily show that 



where, again, no summation over index a is intended. Thus the coefficients of the normal curvature Q.a = 
^a/Sai^x^ are given by ^ai3a = Jj^"" c^^hap. It is clear that ^ajSa = if and only if /,„ ~ 0. Therefore, in 
view of the above we conclude the following: timelike or spacelike geodesies can be confined to the submanifold 
if and only if the normal curvature of vanishes. It is interesting to ask what happens in the case of 
null geodesies. Now since for a curve 7 to be a null geodesic of M"^ e needs to be zero, equations (7) and (8) 
imply that 7 would also be a geodesic of M, irrespective of the value of /,a, and hence irrespective of the 
value of the normal curvatures of M*. Thus the motion of photons in this setting is not constrained by the 
extrinsic curvature of the hypersurface. 

The discussion of the geodesies in the previous section allows this theorem to be generalised to cases 
where the geometry of the ambient space M is Riemannian and the hypersurface under consideration is 
replaced by a submanifold of codimension n. This can be readily seen from Eqs. (8) by noting that the 
condition for the vanishing of the extrinsic (normal) curvatures in this case is given by f^a = which ensures 
that the 4D part of the geodesic Eqs. (8) is geodesic. 

The above results give the conditions for the geodesies of the higher dimensional space M to be confined 
to the co-dimension n hypersurface, but do not give any information concerning the stability of such a 
confinement. We shall consider this question in the following sections, using a phase space analysis. This 
generalises a recent result of the authors concerning the embeddings of hypersurfaces with codimension one 
[11]. 



4 Motion in the extra dimensions: an analysis of the phase space 

Defining ^ = allows the equations (8) to be expressed as an autonomous dynamical system given by 



dX 

where P'^{z, y) = — f'"'{i + kbcZ^'z'^) —T'^^z^z'^, with e = 0, 1. Writing the equations as a dynamical system 
allows the equilibrium points (given here by ^ = = ^) as well as their stabilities to be studied. These 
in turn allow a great deal of information to be gained regarding the types of behaviour allowed by the system 
(see e.g. [14]). 

5 Motion of massive particles in the neighbourhood of the sub- 
manifold M^ 



Let us start with the investigation of the motion of massive particles. For simplicity we shall first consider the 
case where the metric kab is Euclidean, i.e. kai, = diag{+l, ... + 1). It is then easy to see that the equilibrium 
points exist if the simultaneous equations /_o = have real roots (which we denote by j/o = {yl, ?/"))• "^^^ 
equilibrium points are then given by (z° = O. y = yo). 

The nature as well as the stability of these equilibrium points can be obtained by linearising equations 
(9) and studying the eigenvalues of the corresponding Jacobian matrix about the equilibrium points. In this 
case this is a 2n x 2n matrix given by 



n = 



0, 

{fab)n 



nxn Inxn 
On, V 77, 



(10) 



where f^ab = gyagyb ) and 0„xn and Inxn are zero and unit nxn matrices respectively. It is not difficult to 
see that the determinant and the trace of O (detfi and I) are given, respectively, by detfi = (— l)"det 
and / = 0. 

For general forms of the warping function the eigenvalues of the system will satisfy a polynomial of order 
2n which would be difficult to analyse analytically. Our primary aim here, however, is sec whether the 
system is in principle capable of providing confinement of particles in the neighbourhood of the spacetime 
(3 + 1) hypersurface in the co-dimension n setting. Thus rather than pursuing the general case we shall ask 
whether there are classes of special warping functions / for which such confinement is possible. 



As an example we shall consider the cases where fij computed at the equilibrium points is zero for all 
i j and real positive numbers for i = j. In such cases the matrix simplifies and the eigenvalues (all of 
which turn out to be pure imaginary) can be readily found: 

X = ±i^/fZ, a = 1....2n, (11) 

where no summation is intended over a. This condition would clearly be satisfied for warping functions of 

the type 

/ = Ec-(yT, (12) 

a 

where n is a positive integer. 

Writing the system (9) with (10) in terms of normal coordinates, allows it to be expressed in terms of 
de-coupled harmonic oscillators: 

^-".V, (13) 

where cOj = fjj and again no summation is intended over j. 

Thus in such co-dimension n cases, one can view the motion in extra dimensions as confined to n-tori. 
This generalises the case of the centre equilibrium point that was found in the co-dimension one case recently 
[11] and amounts to a toroidal confinement of particles in the neighbourhood of the spacetime hypersurface. 



6 Confinement through an effective potential 

An analysis of the motion in the extra dimensions can also be carried out by reducing the problem to the 

motion of a particle subjected to the action of an effective potential V = V{y). For generality in this section 
we shall assume that the metric kab in Eq. (1) is a Riemannian (positive-definite) metric. 
The geodesic equations for the extra dimensions part are given by 

t-it.,.»V + /,.e-.„,^|i=0. (14) 
In these coordinates equation (6) becomes 



2. dx-dx^ dy'^dy^ 

which again allows the decoupling of the motion in the extra dimensions with equations 

^ - \kbc,ay''r + f,a (e + kbcy'y") = 0. (16) 
A first integral of the above equation may be found by multiplying (16) by the factor 2y"e^-^ to give 

2e^^ ^ - ^rfc6c,a/r) + 2re^^/,a (e + kb^y^y^) = 0, 



which, in turn, gives 



Integrating we obtain 



9f d .h.n\ de^^ /, .h.,.\ de^^ 



e^^ {kb£y'^)=C-ee^i. 
Now, if we assume that /(O, 0) = 0, then it follows that 

C = + e, 

where y^ = kbcyoVo related to the initial kinetic energy corresponding to the motion in the extra dimensions. 
Thus, we have 

e^f{kbcy'r)=yl-Viy), (17) 
where V {y) can be treated as an effective potential energy given by 

V = e\e^f - 1). 



Given that the metric kab is positive-definite the motion is not allowed in the region where V{y) > yg. On the 
other hand, the particle may be bound to S in the neighbourhood of y = if y = is a point of minimum of 
V{y). We note, however, that in cases where the embedding of the submanifold has codimension greater 
than one we can have bounded motion without the particle crossing the brane. 

As a simple example we consider the case corresponding to an embedding of with codimension two, 
that is, M = x iV^. Let us assume the geometry of N'^ is given by the line element 

df = u'^{r){dr^ +r^d9^). (18) 

Clearly wc; can interpret the radial coordinate r as related to the distance between the points of M and 
the submanifold M"*, while the function u{r) gives a measure of the non-Euclideanicity of the metric k^b. 
Also since two-dimensional Riemannian manifolds are conformally flat, any two-dimensional metric kab can 
locally be put in the form (18), but in general with u = u{r, 9) [15]. We now make the assumption, which in 
this case seems to be rather natural, that the warping function / is also a function of r, that is, / = /(r). 
With these assumptions the equation of the motion (17) yields 

e^V {f+rH-") =yl-V{r). 

Prom the equation of motion (16) for 6, we get the following constant of the motion: 

v^r^O = L = const. 
In this way, the equation for the radial motion reduces to 

Here L can be thought of as being related to the angular momentum of the motion of the particle around 
the submanifold M*. Clearly the existence of bounded states depends upon the behaviour of the effective 
potential 



Veff=(V{r)+L' 



We shall end this section by deducing some general properties of the motion. Let us assume that / (0) = 
and that u(0) is regular and docs not vanish at the origin, since otherwise the metric g would not be well 
defined at the spacetime M^. Now if L 7^ then the dominant term in the potential Veff in the limit r ^ 



is L'^^^ — > 00 . We therefore have our first result (which is physically obvious) thus: if i 7^ 0, then the 

particle cannot return to the submanifold, whether or not its motion corresponds to a bound state. On the 
other hand, for L = 0, 14// acts as confining potential (for any value of ^q) if asymptotically (as r — » 00), 
we have / (r) — > 00. 

As a simple application of this result, let us consider a particular case found in the literature [16], in 
which the functions /(r) and u{r) are explicitly given by 

„2/_£!±^ 



e = 



u 



2 



c 



4 



(c2 + r2)2 ' 

where a and c are real constants. In this case the effective potential is given by 



Veff = 



fd, _|_ J.1 



2(c^ + a!:!)(c^+i!) 



The behaviour of 14// is displayed in the Fig. 1, for two different values of the parameter L. The upper 
curve (with L = 1) corresponds to the case where K;// ^ 00 as r ^ and hence the particle cannot cross 
the submanifold M^, despite the fact that it might be in a bound state. The lower curve (with L = 0), on 
the other hand, corresponds to the case where the particle is confined to M*. 



As can be seen, Ve// will act as a confining potential if a > 1 and L = 0. In such cases, there exists a 
limiting value of the kinetic energy above which the particle will escape to infinity. 
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Figure 1: Figure showing the behaviour of the effective potential Vg// for different vahies of L. For the upper 
curve (with L = 1), Vg// ^ oo as r ^ 0, implying that the particle cannot cross the submanifold TW. For 
the lower curve (with L = 0), on the other hand, the particle is confined to M^. The value of the parameter 
c was taken to be c = 100. 

7 Conclusions 

In this paper we have examined some aspects of the motion of massive particles and photons in a (3 + 1 + n)- 
dimensional warped product spaces. Spaces of this type, where the codimension of the embedding is one or 
two, have received a great deal of attention over the recent years mainly in connection with the so-called 
braneworld scenarios. Our treatment has been geometrical and classical in nature. 

We have derived the conditions under which timelike and spacelike geodesies in the full space M coincide 
with those on the codimension n hypersurface. We also have shown that the motion of photons does not 
depend on the extrinsic curvature. Employing the splitting that naturally occurs in such spaces between 
the motion in the hypersurface and the remaining dimensions, and using plane analysis, further allows the 
stability of such a confinement to bo also studied. Using this approach, we have found a novel form of quasi- 
confinement (namely toroidal confinement) which is neutrally stable. The importance of such confinements 
is that they are due purely to the classical gravitational eff'ects, without requiring the presence of brane type 
confinement mechanisms. 

Finally in connection with our results regarding the conditions that need to be satisfied by a warped 
product space in order to ensure the confinement and stability of geodesic motions, we would like to mention 
a recent work [17], in which the author considers an analogous question. However, although this analysis 
is quite general, the author restricts himself to the case of codimension one embedding. Both approaches 
would essentially lead to same results when the warped space is taken to be five-dimensional. 
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